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Abstract 

The quasiparticle excitations in the mixed state of Cd2Re207 have been 
studied by means of muon spin rotation/relaxation (/zSR). The tempera- 
ture dependence of the magnetic penetration depth (A) is consistent with 
a nearly isotropic superconducting order parameter, although a slight dis- 
crepancy which is dependent on the details in the analysis may be present. 
This is also supported by the relatively weak field dependence of A. 
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A class of metal oxides isostructural to mineral pyrochlore has been attracting consid- 
erable attention because they exhibit a wide variety of interesting physical properties. |I[ 
The pyrochlore has a general formula of A2B2O7, consisting of BOq octahedra and eightfold 
coordinated A cations, where A and B are transition metals and/or rare-earth elements. In 
particular, the B sublattice can be viewed as a three-dimensional network of corner-sharing 
tetrahedra, providing a testing ground for studying the role of geometrical frustration in sys- 
tems which have local spins at B sites with antiferromagnetic (AFM) correlation. || Such 
systems are known to remain frustrated even when the exchange interaction is ferromagnetic 
(FM), provided that the spin correlation has local Ising anisotropy. Recent studies have 
revealed a rich variety of phenomena seemingly related to the geometrical frustration, such 
as the occurrence of a spin-glass (SG) phase in i? 2 Mo 2 7 with R—Y, Tb and Dy @J^], the 
unusual behavior of ordinary and anomalous Hall coefficients in the same compound with 
R=Nd, Sm and Gd and the "spin-ice" phase in i2 2 Ti 2 7 with R=Dy and Ho. ||| 

Although metallic pyrochlores comprise a minority subgroup of the pyrochlore family, 
they consist of distinct members such as Tl 2 Mn 2 07, which exhibits colossal magnetoresis- 
tance. PJTUfl Moreover, the recently revealed superconductivity in Cd 2 Re 2 07, a 5d transition 
metal pyrochlore |]TT| , |r2|| , demonstrates that the pyrochlores provide a fertile field for elec- 
tronic correlation adjacent to the perovskite compounds. In this context, it is noteworthy 
that LiV 2 04, a cubic spinel compound in which the V sublattice is isostructural to the B 



sublattice in pyrochlore, behaves similarly to a heavy fermion metal. |T3] 

It is reported that Cd 2 Re 2 07 falls into the bulk superconducting state below T c ~ 1 ~ 
2 K, as confirmed by a large jump of specific heat AC C as well as large diamagnetism 



due to the Meissner effect associated with the occurrence of zero resistivity. [ 1 1 1 The dc 



magnetization curve indicates that the superconductivity is of type II with the upper critical 



field close to 0.29 T at K. [14]] The ratio lSC c j^jT c (with 7 being the Sommerfeld constant) 
is reported to be 1.15, which is smaller than the predicted value of 1.43 for isotropic BCS 
superconductors. Unfortunately, these measurements were performed only above 0.4 K and 
are thus inconclusive in determining the detailed characteristics of superconductivity in 



Cc^R^CV In this Letter, we report on the quasiparticle excitations in the mixed state of 
Cd 2 Re 2 7 studied by muon spin rotation/relaxation (/xSR). The magnetic penetration depth 
A, which reflects the population of normal electrons ( "quasiparticles" ) in the superconductive 
state, is determined microscopically by measuring the muon spin relaxation due to the 
spatial inhomogeneity of magnetic induction in the flux line lattice (FLL). We show that 
the temperature dependence of A is more or less consistent with the prediction based on the 
BCS superconductors with an isotropic gap, although a slight discrepancy is suggested by 
the detailed analysis. 

The /xSR experiments on both single-crystal and polycrystalline Cd2Re2C>7 were per- 
formed on the M15 beamline at the TRIUMF muon facility which provides a beam of nearly 
100% spin-polarized positive muons of momentum 28.6 MeV/c. The specimen was mounted 
on the coldfinger of a 3 He- 4 He dilution refrigerator and cooled from a temperature above 
T c after setting the magnetic field at every field point (i.e., field-cooling) to eliminate the 
effect of flux pinning. The field and temperature scan data were obtained at T = 0.2 K and 
at H = 0.1 T, respectively. Muons were implanted into the specimen (measuring about 10 
mmxlO mm and 1 mm thick) after being passed through a beam collimator. The initial 
muon spin polarization was perpendicular to the magnetic field H and thus to the FLL in 
the superconducting state. 

Since the muons stop randomly along the length scale of the FLL, the muon spin pre- 
cession signal P(t) provides a random sampling of the internal field distribution B(r), 



where 7 M is the muon gyromagnetic ratio (=2nx 135.53 MHz/T), and n(B) is the spectral 
density for muon precession determined by the local field distribution. These equations indi- 
cate that the real amplitude of the Fourier-transformed muon precession signal corresponds 
to the local field distribution n(B). The local field distribution can be approximated as the 
sum of magnetic induction from isolated vortices in the London model to yield 




(0.1) 




dB' 



(0.2) 
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where K is a translation of the vortex reciprocal lattice, Bo (~ H) is the average internal 



field, A is the London penetration depth, and £ v is the cutoff parameter. The term 0(K 2 , K 2 ) 
denotes the nonlocal effect in which the electromagnetic response kernel Q(K) generating 
the supercurrent around the vortex depends on K. While this term is eliminated in the 
conventional BCS superconductors with isotropic s-wave pairing, it becomes important for 
the moe complex order parameters such as anisotropic s-wave or <i-wave (e.g., d x 2_ y 2). The 
London penetration depth in the FLL state is related to the second moment (AB 2 ) = 
((B(r) — -Bo) 2 ) of the field distribution reflected in the /^SR line shape (where ( ) means the 
spatial average). In polycrystalline samples, a Gaussian distribution of local fields is a good 
approximation, where 

P(t) ~ exp(-(T 2 t 2 /2) exptiipHt), (0.4) 



a = 1 ^(AB 2 }. (0.5) 

For the case of ideal triangular FLL with isotropic effective carrier mass m* and a cutoff 
K m 1.4/£u provided by the numerical solution of the Ginsburg-Landau theory, the London 
penetration depth A (with 0(K 2 ,K 2 ) = 0) can be deduced from a using the following 
relation 

a [fis- 1 ] = 4.83 x 10 4 (1 - h)\- 2 [nm], (0.6) 
where h = H/H c2 . While the above form is valid for h < 0.25 or h > 0.7, a more useful 



approximation valid for an arbitrary field is [17 



a = 4.83 x 10 4 (1 -h)[l + 3.9(1 - h) 2 ] 1/2 X' 2 [nm]. (0.7) 

In both cases, A is related to the superconducting carrier density n s as 



indicating that A is enhanced upon the reduction of n s due to the quasiparticle excitations. 
For simplicity, we adopt eq. (p.6|) for the following analysis. 

In a preliminary analysis, we found that the spin relaxation rate due to the FLL formation 
is less than 0.1 fis~ l which is typically of the same order of magnitude as that due to 
static random local fields from nuclear magnetic moments. This means that the additional 
relaxation due to m,113 Cd and 185,187 Re nuclear moments must be considered for the proper 
estimation of A in Cd2Re207. To this end, the following equation was used in the actual 
fitting analysis of the time spectra, 

P(t) = exp[-(a t) v + <J 2 t 2 /2)\ exp(il»Ht), (0.9) 

where o"o is the relaxation rate due to the nuclear moments and v is the power of relaxation. 
The parameters o"o and v were evaluated by fitting the time spectra above T c with a set to 
zero, yielding <jq ~ 0.057(1) /zs _1 and v ~ 1.0 in the polycrystalline specimen at H = 0.1 T. 
Then, a due to the formation of FLL below T c was deduced by analyzing data with cr and 
v being fixed to the above values. A similar analysis was performed for the data obtained 
for the single crystals. 

Figure |I|(a) shows the temperature dependence of a in Cd 2 Re 2 07 at H = 0.1 T. Upon the 
onset of FLL formation, a exhibits a gradual increase with decreasing temperature just below 
T c (0.1 T) ~ 0.7 K at this field. According to the empirical two-fluid model approximately 
valid for conventional BCS superconductors, we have 

A(t)=A(0)- 7 =L=, (0.10) 

which leads to 

a(t) = ( x(0)(l-t 4 ), (0.11) 
where t = T/T c (0.1 T). The fitting analysis by the same formula with an arbitrary power, 

cr(t) = ct(0)(1 — t^), (0.12) 
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with T c as a free parameter yields (3 = 2.8(5) and T c = 0.65(3) K. We also found that 
T c = 0.62(1) K when f3 = 4 is assumed, which is slightly lower than the value T c ~ 0.7 K 
estimated from the specific heat measurement for this field. Although the difference is not 
obvious between these cases (solid curve for f3 = 2.8 and the dashed curve for /3 — 4) in 
Fig. 1(a), we note that the reduced x 2 f° r the former is almost two times smaller (better) than 
the latter. Taking H c2 = 0.29 T, the penetration depth extrapolated to T = (A(0, 0.1 T)) 
is 700(8) nm, leading to the Ginsburg-Landau parameter k ~ 21 with £ ~ 34 nm estimated 
from H c2 (0). Thus, it is concluded that Cd2Re207 is a typical type II superconductor with 
large A, which is consistent with the results of magnetization measurement. fTlH 

The finding that (3 ~ 2.8 may indicate that the deviation AA = A(t) — A(0) oc exhibits 
a tendency predicted for the case of line nodes (o?-wave pairing) with some disorder (i.e., 
dirty limit), where A A oc T 2 . [0 Such a temperature dependence has actually been observed 



in high-T c (YBCO) cuprates. [0,0] Compared with the case of isotropic gap A^ = A , the 



quasiparticle excitations are enhanced along nodes (|A^| = 0) to reduce average n s , leading 
to the enhancement of A. However, the result of the fitting analysis also suggests that the 
discrepancy may be attributed to experimental uncertainty including the precise value of 
T c . Thus, we are led to conclude that the order parameter in Cd2Re207 is mostly isotropic 
with a possibility of residual weak anisotropy as suggested by the slightly reduced value of 
(3. 

As shown in Fig. 2(a), a decreases with increasing external field, where the general 
field dependence is determined by the increasing contribution of normal vortex cores and 
the stronger overlap of field distribution around the cores which are described by the term 
(1 — h) in eq.( |0.6|) . The fitting analysis by eq. ( |0.6j ) with H c2 and X(H = 0) as free parameters 
yields H c 2 = 0.37(5) T and \{H = 0) = 796(12) nm at 0.2 K. Since the deduced value of H c2 
is consistent with that obtained from the magnetization measurement (~ 0.29 T), we can 
conclude that the observed field dependence of A is mostly due to the vortex core/overlap 
effect, except for the fields below ~0.06 T where a slightly steeper reduction of a is suggested. 
The penetration depth deduced for each field using eq. (p.6|) with H c2 = 0.37 T is plotted in 
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Fig. 2(b). In order to evaluate the relative strength of the pair-breaking effect, it is useful 
to introduce a dimensionless parameter r\ to describe the field dependence of A with the 
following simple linear relation, 

A = A(0.2K, h)[l + r]h], (0.13) 

where h = H/H c2 (0.2K) with F c2 (0.2K) = 0.37 T. From the analysis of data in Fig. 2(b) by 
eq. ( |0~l3l) , we find that r] = 0.38(14) with A(0.2K, 0) = 741(5) nm for < H < 0.06 T. 

In general, the field dependence of A is enhanced by two different mechanisms, i.e., the 
nonlinear effect in the semiclassical Doppler shift of the quasiparticle energy levels due to 
the supercurrent around the vortex cores, and the nonlocal effect which further modifies 
the quasiparticle excitation spectrum in the momentum space. In particular, the nonlocal 
effect is important in the system with line nodes because the coherence length is inversely 
proportional to the order parameter, such that (,o(k) = hvp/TrA^. The divergence of £o along 
the nodal directions |A^| = means that the response of quasiparticles near the nodes is 
highly nonlocal. 

It is predicted that rj <C 1 for the isotropic s-wave pairing because the finite gap prevents 
the shifted levels of quasiparticle excitations from being occupied at low temperatures. This 
is supported, for example, by the recent observation in CeRu2, in which i] ~ over the 
field region < h < 0.5 where the system behaves more or less as a conventional BCS 



superconductor with isotropic s-wave pairing. J21J On the other hand, stronger field depen- 
dence is predicted for the case of c?-wave or anisotropic s-wave pairing due to the excess 
population of quasiparticles in the region where |A^| is small or zero. Typical examples 
for the <i-wave pairing are those of high-T c cuprates in which 7] is reported to be 5~6.6 for 



YBCO. |^2[ Meanwhile, in the case of an s-wave superconductor YNi 2 B 2 C in which strong 
anisotropy for is suggested experimentally |23| , r\ ~ 1 is reported from a /iSR study. |24j] 
The comparison of these earlier results with our result suggests that the anisotropy of the 
order parameter in Cd2Re2C>7 is considerably smaller than YM2B2C. 

Finally, we comment that the recent observation of a clear coherence peak below T c in 



the 187 Re NQR measurement [25] does not necessarily indicate the absence of anisotropy 
in the order parameter. In the case of the anisotropic energy gap, the magnitude of the 
coherence peak depends on the mean free path I. The coherence peak is enhanced in a 
certain condition of I where quasiparticles probe only a limited region of the Fermi surface. 
Thus, while the actual value of / is difficult to estimate because of the semimetallic character 
of this compound, the NQR result does not completely rule out the presence of anisotropy 
in general. 

In summary, we have investigated the quasiparticle excitations in the mixed state of 
Cd2Re2C>7 by fiSK. The temperature and field dependence of the London penetration depth 
indicates that the basic feature is consistent with the isotropic order parameter for BCS 
s-wave pairing, although there remains a certain subtlety suggested by the small deviation 
from the theoretical prediction which may be better understood by considering a weak 
anisotropy. 
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FIGURES 

FIG. 1. Temperature dependence of a) the muon spin relaxation rate a due to flux line lattice 
and b) the magnetic penetration depth A, where solid circles are the data from the polycrystalline 
specimen and open triangles are those from single crystals. Solid curves are results of fitting by a 
relation a oc 1/A oc 1 — (T/T C )P with (3 and T c being free parameters. Dashed curves are obtained 
when (3 = 4 with T c clS cl free parameter. 

FIG. 2. Magnetic field dependence of a) the muon spin relaxation rate a due to flux line 
lattice and b) the magnetic penetration depth A. The solid line in a) is a result of fitting by 
eq.( p.6| ) with H C 2 being a free parameter. Dashed lines are obtained by fitting data with a relation 
A oc 1 + r](H/H c2 ) only for < H < 0.06 T. 
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